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We explore the intensity correlations for thermal photons having Kt < 2 GeV/c, for central 
collisions of heavy nuclei at RHIC and LHC energies. These photons get competing contributions 
from the quark and the hadronic phases. This competition gives rise to a unique structure, especially 
in the outward correlation function, due to the interference between the photons from the two 
sources. The temporal separation of the two sources provides the life time of the system and their 
strengths provide the relative contribution of the two phases. The results are found to be quite 
sensitive to the quark-hadron phase transition temperature and the formation time of the plasma. 



I. INTRODUCTION 

The last two decades have witnessed a concerted and 
a well co-ordinated theoretical and experimental effort 
to produce and study quark-gluon plasma - the decon- 
fined strongly interacting matter, in relativistic colli- 
sions of heavy nuclei. The eminent commissioning of 
the Large Hadron Collider (LHC) at CERN and various 
upgrades, both in the accelerator and the detection sys- 
tems, at the Relativistic Heavy Ion Collider (RHIC) at 
Brookhaven will offer further opportunities to advance 
our understanding of the physics of this novel state of 
matter, which permeated the early universe. We have 
already been rewarded with the discovery of the jet- 
quenching [l|, 13] , the elliptic flow 0, 0| , and the successful 
formulation of the partonic recombination Q as a model 
for hadronization in recent experiments at RHIC. The 
thermal photons radiated from such collisions are being 
studied to determine the high temperatures reached in 
these collisions [1, 0, E H • 

The next step in this endeavour involves a critical ex- 
amination of our concepts about the formation and the 
evolution of the plasma. An important question in this 
connection is how quickly, if at all, does the plasma ther- 
malize; so that we can use the powerful methods of hydro- 
dynamics to model the evolution of the system. Theoret- 
ical estimates for the thermalization time (to) vary con- 
siderably and experimental observables like elliptic flow 
of hadrons only suggest that tq ~ 1.0 fm/c. A more quan- 
titative experimental determination of To would be very 
valuable. One would also like to know the lifetime of the 
interacting system. Can we determine the temperature 
at which the phase transition takes place? 

The quantum statistical interference between identi- 
cal particles emitted from these collisions is expected to 
provide valuable inputs for the space-time description of 
the system. The use of photons for these studies admits 
several advantages. Firstly, they interact only weakly 
with the system after their production. Thus, they are 
not subjected to distorting effects such as re-scattering 
and Coulomb interactions which affect the results for 
hadrons. Secondly, and even more importantly, photons 
are emitted from every stage of the collision dynamics. 



These aspects give us a hope of getting a direct infor- 
mation about the earliest, hot and dense stage of the 
system by studying the photons having larger transverse 
momenta, Kt- Recently, some additional sources of large 
Kt photons have also been proposed [TH EH • 

The difficulty, of-course, arises from a meager emission 
of direct photons which lie buried in the huge background 
of decay photons. So far only one measurement involv- 
ing photons having very low Kt by the WA98 Collabora- 
tion [12[ for the central collision of lead nuclei at CERN 



SPS has been reported. 

On the other hand, the theory of the intensity interfer- 
ometry of photons from relativistic heavy ion collisions 
has been pursued in considerable detail by several au- 
thors [H, [3, HU . It is generally felt that the experimen- 
tal efforts for these studies have a larger likelihood of 
success at RHIC and LHC energies because of larger ini- 
tial temperature of the plasma and a large suppression 
of pions due to jet-quenching. Of late, there have also 
been tremendous advances in methods for identification 
of single photons [8|. 

In the present work, we focus our attention on the pho- 
tons having intermediate Kt ~ 0.2 - 2 GeV/c. The pho- 
tons having Kt <C 2 GeV/c will mostly originate from 
the hadronic phase of the system. They are expected to 
reveal a source which should be strongly affected by the 
flow and expansion of the system. The photons having 
Kt ^S> 2 GeV/c should unveil a source which is in the 
infancy of the hot and dense quark-gluon plasma, and 
where the flow has just started to develop. The photons 
having Kt < 2 GeV/c are unique. They have their origin 
either in the hot and dense quark phase of the system or 
in the relatively cooler but rapidly expanding hadronic 
phase, where a large build-up of the radial flow boosts 
their transverse momenta. 

We shall see that this leads to a rich structure in the 
correlation function for thermal photons especially when 
studied as a function of the outward momentum differ- 
ence (q , see later), due to the interference of the two 
sources. The two sources also manifest in the correla- 
tion functions for the longitudinal (qi) and sideward (q s ) 
momentum differences. This renders the correlation very 
sensitive to the formation time of the plasma. An early 
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formation and thermalization would provide a large ini- 
tial temperature, while a late formation and thermaliza- 
tion will lead to a smaller initial temperature. This anal- 
ysis can also help to find the fractional contributions of 
the quark and the hadronic phases to the single photon 
spectrum. 

Since the interference mentioned above is controlled 
by the relative contributions from the quark and the 
hadronic phases, which in turn are decided by the critical 
temperature used in the model, we find a unique sensitiv- 
ity of the results to the temperature at which the phase 
transition takes place. 

We discuss the basic formalism in the next section. 
The results for RHIC and LHC energies are discussed in 
sect. Ill and IV respectively. In sect. V we discuss the 
sensitivity of our results to the initial formation time of 
the plasma and the transition temperature. Finally, we 
summarize our findings in sect. VI. 



II. FORMULATION 

One can define the spin averaged intensity correlation 
between two photons with momenta ki and k 2 , emitted 
from a completely chaotic source, as: 



we can write UA 



C(q,K) = l 



1 |J d i xS{x,~K)e 



ix-q 



(i) 



2fd i xS(x,k 1 ) JtfxSix,^) 

where S(x, 1 K) is the space-time emission function, and 

q = k!-k 2 , K= (ki+k 2 )/2 . (2) 

We shall use hydrodynamics to model the evolution of the 
system. The space-time emission function S is approxi- 
mated as the rate of production of photons , EdN / d 4 xd 3 k, 
from the quark and the hadronic phases. 

The interference between the thermal photons from 
quark and the hadronic phases is best studied by writ- 
ing the source function S as Sq + Sh in the numerator, 
where Q and H stand for the two phases respectively, 
and then including either one, or the other, or both. 

We shall discuss the results for the correlation function 
C(q, K) in terms of the outward, sideward, and longitu- 
dinal momentum differences, q Ql q s and q#. Thus writing 
the 4-momentum as of the ith photon, we have 



(k lT cosh 2/i,ki) 



with 



ki = (k iT cosipi,kiTsmipi,kiT sinh 



(3) 



(4) 



where kx is the transverse momentum, y is the rapidity, 
and is the azimuthal angle. Defining the difference and 
the average of the transverse momenta, 

q T = kiT -k 2 T ,K T = (kiT + k 2T )/2, (5) 
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The radii corresponding to the above momentum dif- 
ferences are often obtained by approximating the corre- 
lation function as, 

C(q ,q s , q t ) = 1 + 1 exp [- (q 2 R 2 + q 2 R 2 + q 2 R 2 )] . (9) 

We also define the root mean square momentum differ- 
ence (qf) as, 



J(C-l)qfd qi 

f(C-l)d qi ' 



(10) 



so that for the Gaussian parameterization given in Eq. [9l 
we have 



R 



1 



2<«?> 



(11) 



Thus, l/[2 {qf)} 1 / 2 becomes a useful measure when the 
correlation function has a more complex nature as we 
shall see later. 

We consider central collision of gold and lead nuclei, 
corresponding to the conditions realized at the RHIC 
and LHC, respectively. We assume that a thermally and 
chemically equilibrated quark-gluon plasma is produced 
at an initial time tq. We further assume an iscntropic ex- 
pansion of the system to estimate the initial temperature, 
Tq, in terms of particle rapidity density. Thus, 



2tt 4 



1 dN 



45C(3) A T dy 



= 4aT J r , 



(12) 



where At is the transverse area of the system, dN/dy 
is the particle rapidity density, and a = 42.257r 2 /90 for 
a plasma of massless quarks (u, d, and s) and gluons. 
The number of flavours for this purpose is taken as w 2.5 
to account for the mass of the strange quarks. For the 
present work we shall mostly consider tq — 0.2 fm/c and 
give illustrative results for tq varying from 0.2 to 1.0 fm/c, 
keeping dN/dy (or the total entropy) fixed. We add that, 
for study of thermal photons a value close to 0.2 fm/c 
may be more appropriate. The initial energy density is 
taken as a weighted sum of wounded nucleon and bi- 
nary collision distributions as in earlier studies [16]. The 
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FIG. 1: (Colour on-line) (a) The outward, (b) sideward, 
and (c) longitudinal correlation functions for thermal photons 
produced in central collision of gold nuclei at RHIC taking 
To = 0.2 fm/c. Symbols denote the results of the calculation, 
while the curves denote the fits. 



quark-hadron phase transition is assumed to take place 
at a temperature of 180 MeV, while the freeze-out takes 
place at 100 MeV. The relevant hydrodynamic equations 
are solved under the assumption of boost-invariant longi- 
tudinal and azimuthally symmetric transverse expansion 
using the procedure discussed earlier [l~7j and integration 
performed over the history of evolution. We use the com- 
plete leading order results for the production of photons 
from the quark matter [l8j], and the results from Tur- 
bide et al. [191 ] for radiation of photons from the hadronic 
matter. A rich equation of state with the inclusion of all 
the particles in the particle data book, having M < 2.5 
GeV/c 2 describes the hadronic matter. We take dN/dy 
at y = as 1260 [nlfor 200 A GeV Au+Au collisions 
at RHIC and 5625 ^20j for 5.5A TeV Pb+Pb collisions 
at LHC. A smaller value for the particle rapidity den- 
sity at LHC may perhaps be more appropriate, though. 
This will, however, not change the nature of the findings 
reported here. 



III. RESULTS FOR RHIC ENERGIES 

As a first step, in Fig. [T] we have shown the outward, 
sideward, and longitudinal correlation functions for ther- 
mal photons at RHIC having K T « 2 GeV/c. The 4- 
momenta of the two photons are chosen so that when we 
study the outward correlations, the q s and qe are identi- 
cally zero and the dependence on q is clearly seen, and 
so on. We first discus the results when only the hadronic 
matter contribution or only the quark matter contribu- 
tion is included in numerator (Eq. [T]). We find that the 
correlation functions for the two phases can be approxi- 
mated as: 



C{q h a) = 1 + 0.5|p ija | 5 



(13) 



where i = o, s, and t, and a denotes quark matter (Q) 
and hadronic matter (H) in an obvious notation. The 
source distribution \pi. a \ is very well described by, 



\Pi,a\ 



Ii exp [- 0.5 (q 2 R 2 ,J] 



where, 



and 



dNr 



(dN Q + dN H ) ' 



Ih = 



dN 



H 



(dN Q + dN H ) ' 



(14) 



(15) 



(16) 



The final correlation function denoted by 'sum' in the 
figures can be approximated as: 



Cfe) = 1 



0.5 [ |p^ Q | 2 + Kh| 2 

2 \pi,Q\\pi,H\ cos(qARi) 



(17) 



which clearly brings out the interference between the two 
sources [2l|. Here ARi stands for the separation of the 
two sources in space-time and q is the 4-momentum dif- 
ference. For thermal photons having Kt ~ 2 GeV/c at 
RHIC, various radii (in fm) are obtained as: 



Ro,Q — 2.8, R ,H 
Rs,Q ~ R s ,h = 2i 
Rt,Q = 0.3, Re,H 



7.0, AR 
AR S « 0., 
1.8, AR e i 



12.3, 



0. 



(18) 



These results imply [21( that, while the spatial sep- 
aration of the two sources is negligible, their temporal 
separation is about 12 fm. This gives the life-time of 
the system. If the mixed phase is of shorter duration or 
absent, this will obviously decrease. 

This is seen more clearly in Fig. [5] where we have 
plotted the Kt dependence of AR a and the outward, 
sideward, and longitudinal radii for the hadronic and 
the quark-matter sources of photons, obtained using the 
above procedure. We see that the outward, sideward and 
the longitudinal radii for the quark contribution depend 
weakly on the transverse momentum, indicative of only 
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FIG. 2: (Colour on-line) Transverse momentum dependence 
of (a) the outward radii and temporal duration along with (b) 
the sideward and (c) the longitudinal radii for the hadronic 
and quark-matter sources obtained by fitting the final cor- 
relation function for thermal photons at RHIC energy. The 
radii determined from the root mean square momentum dif- 
ference for the correlation function are also given for a com- 
parison p^ |. 



a mild development of the flow during that phase. The 
corresponding radii for the hadronic contribution show a 
stronger dependence on the transverse momentum, which 
is indicative of a more robust development of the radial 
flow during the hadronic phase. The duration of the 
source reveals a very interesting structure, and in fact 
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FIG. 3: (Colour on-line) Transverse momentum dependence 
of fraction of thermal photons from quark matter (Iq) and 
hadronic matter (Ih) at RHIC energy. The solid curves gives 
the results for Tc = 180 MeV, while the dashed and dot- 
dashed curves give the results for Tc equal to 160 and 200 
MeV respectively, (see text). 



it rises slightly at higher transverse momenta as the pho- 
tons emitted during the hadronic phase benefit from the 
strong radial flow to greatly increase their transverse mo- 
menta. The saturation in AR a towards low Kt has its 
origin in the competition between radial expansion and 
the decoupling of hadronic matter as it cools down below 
the freeze-out temperature at the edges. 

The inverse root mean square momentum which, as we 
commented earlier, is a measure of the correlation radius, 
is seen to vary rapidly with Kt, for the outward correla- 
tion due to rapid variation of the competing contributions 
of the quark matter and hadronic matter phases. For the 
sideward correlation, it goes over smoothly from a value 
which is close to that for the hadronic matter at lower 
Kt to that for the quark matter at higher Kt, as one 
would expect. (The slight difference of l/[2 (q 2 )] 1 ^ 2 from 
R s ,q ~ R s ,h at large Kt is due to the deviation of the 
calculated correlation function from a perfect Gaussian, 
to which it is fitted.) 

As the fractions of the quark matter and the hadronic 
matter contributions play such a pivotal role in the mo- 
mentum dependence of overall correlation function, wc 
study them next. The fractions Iq and Ih taken from 
the calculations for the single photon spectra are shown 
in Fig. [31 along with their dependence on the critical 
temperature. We shall discuss the importance of this 
dependence a little later. 

In order to understand these interesting results more 
clearly, we have plotted the source- functions as a function 
of the transverse distance and time for thermal photons 
having K T « 0.5, 1.0, and 2.0 GeV/c, in Fig.g]and Fig.[5] 
In all the cases we see that the radial distributions for the 
QGP as well as for the hadronic phase are centered at tt 
near zero and that the source function for the hadronic 
phase extends considerably beyond the same for the QGP 
phase. This of-course is only to be expected, due to the 
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FIG. 4: (Colour on-line) The radial dependence of the source- 
distribution function for emission of photons having trans- 
verse momenta of (a) 0.5, (b) 1.0, and (c) 2.0 GeV/c at RHIC 
energy. 



FIG. 5: (Colour on-line) The temporal dependence of the 
source distribution function for radiation of photons having 
transverse momenta of (a) 0.5, (b) 1.0, and (c) 2.0 GeV/c at 
RHIC energy. 



large transverse expansion of the system. We also note a 
good yield of thermal photons from the hadronic phase 
at larger radii. In fact, one can notice a slightly enhanced 
emission at larger radii compared to that from the cen- 
tral region for Kt — 2 GeV/c. This, as we noted earlier, 
arises due to the large kick received by the photons due 
to the radial expansion, leading to the blue-shift of their 
transverse momenta. The relative importance of the two 
contributions will depend on Kt and the transition tem- 
perature and lead to a rich structure in the resulting cor- 



relation function. However, note that the second term in 
the source of photons from the hadronic matter becomes 
important at larger Kt, where the contribution of the 
hadronic phase to photons is rather small, and thus it 
could be difficult to detect this effect. 

Even a more valuable insight is provided by the tem- 
poral structure (Fig. [5]) of the emission of photons from 
the QGP and the hadronic sources; the latter emerging 
after a lapse of some time, which the system spends in the 
QGP phase. The duration over which the photons from 
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FIG. 6: (Colour on-line) Same as Fig. Q] for Pb+Pb at LHC. 



the hadronic phase are emitted is quite large and essen- 
tially controls the parameter AR Q . It is worth recalling 
that the temporal structure of the source function seen 
here is qualitatively similar to the one seen in nucleus- 
nucleus collisions at cyclotron energies plj . 



IV. RESULTS FOR LHC ENERGIES 
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The Large Hadron Collider will study collisions of lead 
nuclei at unprecedented energies of 5500 GeV/A. It is ex- 
pected that the initial temperature likely to be attained 
reached in such collisions would be much higher than 
that at RHIC energies. This would provide a golden op- 
portunity to study the properties of the QGP and the 
dynamics of its evolution. A larger initial temperature 
would lead to a larger duration of the interacting system, 
which in turn will provide an ample opportunity to the 
mechanism of expansion to develop. We can thus expect 
to put our models about the evolution of the interacting 
system to a very rigorous test. 

In Fig. [6] we show our results for the intensity cor- 
relation of thermal photons having Kt ~ 1-5 GeV/c 
produced at LHC. We see, as before, an interference of 
the photons from the hadronic and QGP phases of the 



FIG. 7: (Colour on-line) Transverse momentum dependence 
of (a) the outward radii and temporal duration along with (b) 
the sideward, and (c) the longitudinal radii for the hadronic 
and quark-matter sources obtained by fitting the final corre- 
lation function for thermal photons at LHC energy. The radii 
determined from the root mean square momentum difference 
for the correlation function is also given for a comparison [l3j]. 



system. We add that the fit to the calculated values 
for the outward correlation function can be improved 
considerably (see dashed curve) by adding one more 
Gaussian term, centered at go ~ 0.06 GeV/c in the 
source term for the hadronic matter, or by approximat- 
ing p H ~ \/^{Ch — 1), where Ch is the corresponding 
correlation function obtained numerically. We are trying 
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FIG. 8: (Colour on-line) Transverse momentum dependence 
of fraction of thermal photons from quark matter (Iq) and 
hadronic matter (Ih) at LHC energy. The solid curves gives 
the results for Tc =180 MeV, while the dashed and dot- 
dashed curves give the results for Tc equal to 160 and 200 
MeV respectively. 



to understand this observation. 

The results for C(q s ) and C(qe) are similar in nature 
to those found at RHIC energy. We note that now the 
temporal separation of the two sources is about 19 fm/c 
(see inset Fig. [5]). 

Next we discuss the Kt dependence of the correla- 
tion radii (Fig. [7]). We see a behaviour which is qualita- 
tively similar to the one obtained earlier, though all the 
correlation radii are larger, specially AR which gives 
the duration of the source. We do emphasize here that 
l/[2 (g 2 )] 1 / 2 is closer to R^q at LHC due to the domi- 
nance of quark matter contribution even at modest Kt- 

In Fig [51 we have plotted the fractions of momen- 
tum dependence of quark and hadron contribution at 
LHC energies, at three transition temperature Tc {= 
160, 180, and 200 MeV). Comparing these results with 
those of Fig. [3j we see once again that decreasing Tc in- 
creases the fraction of photons coming from the quark 
matter. We also note that the transverse momentum 
where the quark and hadronic contributions become 
equal is quite sensitive to the transition temperature Tc- 
Realizing that, in an ideal situation we would be able to 
decompose the outward, sideward, and longitudinal cor- 
relations into two sources and that their intercepts on the 
y-axis will give Iq and Ih, this opens up the tantalizing 
possibility of determination of the transition tempera- 
ture, the two fractions becoming equal at Kt- 

We next discuss the spatial (Fig. O and temporal 
(Fig. [T0f distribution of source function for photons hav- 
ing Kt = 0.5, 1.0, and 2.0 GeV/c produced in central 
collision of Pb nuclei at LHC energy. Even though we 
note that the source distributions at LHC are qualita- 
tively similar to those at RHIC, the extension of the 
hadronic sources beyond those for the quark sources 
(Fig. [5]) and enhanced production at larger vt from 
hadronic sources are amplified considerably here. The 
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FIG. 9: (Colour on-line) The radial dependence of the source- 
distribution function for emission of photons having trans- 
verse momenta (a) 0.5, (b) 1.0, and (c) 2.0 GeV/c at LHC 
energiy. 



dying out of the quark source and a delayed build-up of 
hadronic source in time (Fig. llQ[) are seen to emerge very 
clearly. 



V. SENSITIVITY TO r AND T c 

We have seen that the final correlation is decided by 
the relative contributions from the quark-matter which 
occupies a smaller volume and is shorter lived, and those 
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FIG. 10: (Colour on-line) The temporal dependence of the 
source distribution function for radiation of photons having 
transverse momenta of (a) 0.5, (b) 1.0, and (c) 2.0 GeV/c at 
LHC energy. 



from the hadronic matter, which occupies to a larger vol- 
ume and lives longer. 

If the system thermalizes quickly, the initial tempera- 
ture would be large. One can test the sensitivity of the 
results to the formation time of the plasma, by consid- 
ering systems with identical entropy, but having varying 
formation times (tq). The fractional contribution of the 
quark matter, (Iq) will increase with decreasing tq. In 
Fig. 1111 we have shown the To dependence of the outward 
and longitudinal correlations at RHIC. The sideward cor- 
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FIG. 11: (Colour on-line) r dependence of (a) the outward 
and (b) the longitudinal correlation function at RHIC. 
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FIG. 12: (Colour on-line) Tc dependence of the outward cor- 
relation function at RHIC. 



relation function did not show any perceptible change due 
to the variation of To and is not shown here. We must 
add here that choosing a large To would necessitate the 
inclusion of the pre-equilibrium contribution to the pho- 
tons which must surely be there, at least at larger 
Kt- The jet-conversion mechanism [Tl| is also likely to 
contribute at larger Kt- These contributions, per- force 
have their origin in the deconfined matter and their spa- 
tial distribution is not likely to be very different from 
the early stages of the initial distributions assumed here. 
Their contributions would increase Iq in these studies, 
which in turn could mimic an effectively smaller tq . This 
would still be useful, as it will amount to getting an ef- 
fective tq after which the hydrodynamics can be applied. 
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However these discussions open the door to another in- 
teresting and potentially powerful observation, perhaps, 
with a far reaching implication. We have already noted 
the sensitivity of the fractions of the contributions of 
the quark-matter (Iq) and the hadronic matter (Ih) to 
the transition temperature at both RHIC and LHC (see 
Figs. [3] and [5]). In Fig.[T2]we show our results for the sen- 
sitivity of the outward correlation at RHIC to the transi- 
tion temperature. Recalling that an increase in the tran- 
sition temperature leads to a decrease in Iq, and that the 
quark matter contribution has a smaller R a the change 
in the interference pattern seen at larger q Q is easily un- 
derstood. 

It is felt that the results shown here should be valid 
whenever the correlations arise from contributions from 
two sources separated in space and time (see e.g., 
Ref. [2l|). The contribution of the two sources to the 
correlation function also provides a natural explanation 
for the failure of the earlier studies [13( to find a simple 
Gaussian parametrization for it. 



interferometry of thermal photons at RHIC and LHC en- 
ergies is calculated. The correlation functions are marked 
by a very distinctive interference between the photons 
from the quark and the hadronic matter, which is most 
clearly visible in the outward correlation. We have calcu- 
lated the transverse momentum dependence of the cor- 
relation radii and the duration of the emission and tried 
to understand their behaviour by calculating the spatial 
and the temporal distribution of the sources and their 
contributions. The study has thrown open an interest- 
ing possibility of determination of transition temperature 
and the formation time of the plasma. Finally, we would 
like to add that even though several earlier studies have 
also talked of the two sources of photons (quark and 
hadronic matter) contributing to the correlation func- 
tion, the present work, as far as we know, is the first 
attempt to study their interference in relativistic heavy 
ion collisions. 



VI. SUMMARY 
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